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I. INTRODUCTION 

Since the 1960s the europium chalcogenides EuX 
(X=0, S, Se, Te) have attracted tremendous research ac- 
tivity, experimentally as well as theoretically.Emu They 
are magnetic semiconductors, which crystallize in the 
rocksalt structure with increasing lattice constants (5 to 
7 A) when going from the oxide to the telluride. The Eu 2+ 
ions occupy lattice sites of an fee structure so that each 
ion has twelve nearest and six next-nearest Eu-neighbors. 

As to their purely magnetic properties the EuX are 
considered almost ideal realizations of the Heisenberg 
model for the so-called local-moment magnetism. Their 
magnetism is due to the half-filled 4f shell of the Eu 2+ . 
The 4f charge density distribution is nearly completely 
located within the filled 5s 2 and 5p e shells so that the 
overlap of 4f wavefunctions of adjacent Eu 2+ ions is neg- 
ligibly small. Hunds rules of atomic physics can be ap- 
plied yielding an 8 S7/2 ground state configuration of the 
4f shell. The 7fj,B moments are exchange coupled re- 
sulting in antiferromagnetic (EuTe, EuSe), ferrimagnetic 
(EuSe), and ferromagnetic (EuO, EuS) orderings at low 
temperatures. The fact that the magnetic contribution to 
the thermodynamics of the EuX is excellently described 
by the Heisenberg Hamiltonian 
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allows to test models of the microscopic coupling mecha- 
nism by direct comparison to experimental data. There 
is convincing evidence that the exchange integrals can 
be restricted to nearest (Ji) and next nearest neigh- 
bors (J2)™- Ji is positive (ferromagnetic) decreasing 
from the oxide to the telluride. Ji is negative, except 
for EuO, where the antiferromagnetic coupling increases 
in magnitiide.|£rom the sulfide to the telluride. Liu and 
coworkersEoElEI have proposed an indirect exchange be- 
tween the localized 4f moments mediated by virtual ex- 
citations of chalcogenide-valence band (p) electrons into 
the empty Eu 2+ (5d) conduction bands together with a 
subsequent interband exchange interaction of the d elec- 
tron (p hole) with the localized 4f electrons. Using 



this picture— very similar to the Bloembergen- Rowland 
mechanisrr£3, the calculated values for J\ , J2 agree nicely 
with experimental data, both in sign as well as in magni- 
tude, for EuO, EuS, and EuSe. The results are found by 
a perturbative calculation of the indirect 4f-4f exchange 
interaction (|l|) with data from a linear combination of 
atomic orbitals (LCAO) method as band structure in- 
put. The different distances of the 4f moments do obvi- 
ously create the different magnetic behavior of the EuX. 
For the, exchange integrals of the two ferromagnets one 
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Although in EuO the ferromagnetic interaction is more 
pronounced (T c (EuO) = 69.33K; T c (EuS) = 16.57K)tl a 
greater variety of experiments has been carried out with 
EuS than with EuO. The reason is that si ngle crystals as 
well as films with well defined thicknessescJi3 can better 
be prepared for EuS than for EuO. Apart from this, the 
ferromagnetism of EuS is interesting in itself for two rea- 
sons. There are competing exchange integrals J\ and J 2 , 
and the magnitude of the dipolar energy is comparable 
to the exchange energy. 

Besides the purely magnetic properties a striking tem- 
perature dependence of the (empty) conduction bands 
has caused intensive investigation. This was first de- 
tected for the ferromagnetic EuX as red-shift of the opti- 
cal absorption edge (4f-5d) upon cooling below Tqc3- The 
reason is an interband exchange coupling of the excited 
5d electron to the localized 4f electrons that induces the 
temperature dependence of the localized moment system 
into the the empty conduction band states. A further 
striking effect, which is due to the induced temperature 
dependence of the conduction band states, i s a p etal- 
insulator transition observed in Eu-rich EuOEdlEj. The 
Eu richness manifests itself in twofold positively charged 
oxygen vacancies. One of the two Eu 2+ excess electrons, 
which are no longer needed for the binding, is thought 
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to be tightly trapped by the vacancy. Because of the 
Coulomb repulsion, the other electron occupies an im- 
purity level fairly close to the lower band edge. With 
decreasing temperature below T c the band edge crosses 
the impurity level thereby freeing impurity electrons. A 
conductipiy jump as much as 14 orders of magnitude was 
observedll3. Other remarkable effects result from the in- 
teraction of the band electron with collective excitations 
of the moment system. One of these is the creation of 
a characteristic quasiparticle {magnetic polaron) which 
can be identified as a propagating electron dressed by a 
virtual cloud of excited magnons. 

In previous papers we have proposed a method for 
the determination of the ternoerature dependent elec- 
tronic structure of bulk EuOt£l as well as EuO-filmsEZl. 
The treatment is based on a combination of a multi- 
band Kondo-lattice model (MB-KLM) with first princi- 
ples tight-binding linear muffin-tin orbital (TB-LMTO) 
band structure calculations. The many-body treatment 
of the (ferromagnetic) KLM was combined with the first- 
principles part by strictly avoiding a double-counting of 
any relevant interaction. The most striking result con- 
cerned the prediction of a surface state for thin EuO(lOO) 
films, the temperature shift oLwhich may cause a surface 
halfmetal-insulator transitions. For low enough temper- 
atures the shift of the surface state leads to an overlap 
with the occupied localized 4f states. Therefore, one can 
speculate that the resistivity of the EuO(lOO) films might 
be highly magnetic field dependent, so that a colossal 
magnetoresistance effect is to be expected. 

In this paper we investigate in a similar manner the 
other ferromagnet EuS, where we restrict ourselves first 
to the bulk material. We want to derive the temper- 
ature dependent quasiparticle band structure (Q-BS), 
in particular concentrating on those effects, which are 
due to a mutual influence of localized magnetic 4f states 
and itinerant, weakly correlated conduction band stafaaa 
There was earlier work on the Q-BS of bulk EuSE£lEj. 
In these papers, however, an approach was employed 
that decomposes the Eu-5d band into five consecutive 
non-degenerate subbands. For each of the subbands 
a single-band KLM was evaluated therewith disregard- 
ing the full multiband-structure of the EuS conduction 
band. Obviously this procedure leads to an overesti- 
mation of certain correlation effects as a consequence 
of certainly too narrow subbands. We therefore use in 
this paper a multiband 4f-5d Kondo-lattice model to get 
reliably the temperature dependent Q-BS of EuS with 
all correlation effects in a symmetry-conserving manner. 
Our method combines a many-body analysis of the men- 
tioned multiband-model with a self-consistent LMTO 
band structure calculation. 

Since the technical details can be taken fromlll we 
present in the following only the general procedure to- 
gether with those aspects which are vital for the under- 
standing of the new EuS results. In Section 2 we formu- 
late the multiband Hamiltonian and fix its single-particle 
part by a realistic band structure calculation. Further- 



more, we describe the parameter choice for the decisive 
interband exchange coupling. Section 3.1 is devoted to 
the local-moment ferromagnetism of EuS, while Section 
3.2 repeats shortly how we solved the multiband Kondo- 
lattice model. In Section 4 the electronic EuS structure 
is discussed in terms of quasiparticle band structures and 
densities of states (Q-DOS) as well as spectral densities, 
which are closely related to the angle and spin resolved 
(inverse) photoemission experiment. 



II. MULTIBAND KONDO-LATTICE MODEL 

The complete model-Hamiltonian for a real substance 
with multiple conduction bands consists of three parts: 

H = H 5d + H if + H df (4) 

The first term contains the 5d conduction band struc- 
ture of the considered material as, e.g., EuS: 

H$d = ^ ^ Tj™ m c\ mcr Cj m > a (5) 

i,j ra.ra' 

The indices m and m' refer to the respective 5d sub- 
bands, i and j to lattice sites. c\ mr7 and Cj mcr are, re- 
spectively, the creation and annihilation operator for an 
electron with spin a (a =|, J,) from subband m at lattice 
site R,-. Tj™ m are the hopping integrals, which are to 
be obtained from an LDA calculation in order to incor- 
porate in a realistic manner the influences of all those 
interactions which are not directly accounted for by our 
model Hamiltonian. 

Each site is occupied by a localized magnetic mo- 
ment, represented by a spin operator S;. It stems from 
the half-filled 4f shell of the Eu 2+ ion, according to 
Hund's rule being a pure spin moment of 7fJ.B- The ex- 
change coupled localized moments are described by an 
extended Heisenberg Hamiltonian 

ij i 

In the case of EuS the exchange integrals Jjy can be 
restricted to nearest and next nearest neighbors @ . The 
non-negligible dipolar energy in EuS is expressed by a 
single-ion anisotropy Dq. 

The characterizing feature of the normal single-band 
KLM, also called s-f or s-d model, is an intraatomic ex- 
change between conduction electrons and localized spins. 
The form of the respective multiband-Hamiltonian can 
be derived from the general on-site Coulomb interaction 
between electrons of different subbands. It was shown 
irJlaO that in the special case of EuX (half-filled 4f shell, 
empty conduction band) the interband exchange can be 
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written as: 

H df = ~2 J ^2 { S i ( n '™T " n »™i) ( 7 ) 

im 

J is the corresponding coupling constant, and further- 
more: 

n ima = c^ mo .c; mcr ; S = S i iS' 3 ' (8) 

The first term in (Q) represents an Ising type interaction 
while the two others refer to spin exchange processes. 
The latter are responsible for some of the most typical 
KLM properties. 

In order to incorporate in a certain sense all those in- 
teractions, which are not directly covered by the model 
Hamiltonian, we take the hopping integrals from a band 
structure calculation according.ln the TB-LMTO-atomic 
sphere program of AndersenEiJEa. In this method, the 
original Hamiltonian is transformed to a tight-binding 
Hamiltonian containing nearest neighbor correlations, 
only. The transformation is obtained by linearly combin- 
ing the original muffin-tin orbitals to the short ranged 
tight -binding muffin-tin orbitals. The evaluation is re- 
stricted to 5d bands, only. LDA-typical difficulties arise 
with the strongly localized character of the 4f levels. To 
circumvent the problem we considered the 4f electrons as 
core electrons, since our main interest is focused on the 
reaction of the conduction bands on the magnetic state of 
the localized moments. For our purpose the 4f levels ap- 
pear only as localized spins in the sense of H 4f in Eq.([§. 
Figure [l] shows the calculated spin-dependent band struc- 
ture of EuS, where, of course, the 4f levels are missing. 
Clearly, the conduction-band region is dominated by Eu- 
5d states. For our subsequent model calculations it is 
therefore reasonable to restrict the single-particle input 
from the band structure calculations to the Eu-5d part, 
only. The low-energy part in Figure [I] belongs to the 
S-3p states. For comparison we have also performed an 
LDA+U calculation which is able to reproduce the right 
positions of the respective bands. However this method 
suffers from the introduction of adjustable parameters U 
and J being, therefore, no longer a "first principles" the- 
ory. The results of the LDA+U calculation do not differ 
strongly from those of the "normal" LDA, with the 4f 
electrons treated as core electrons. So we have chosen 
the much simpler LDA calculation. Since we are mainly 
interested in overall correlation and temperature effects, 
the extreme details of the bandstructure are surely not 
so important. In Figure ||, from the same calculation, the 
LDA-density of states is displayed. A distinct exchange 
splitting is visible which can be used to fix the interband 
exchange coupling constant J in Eq.(0). Assuming that 
an LDA treatment of ferromagnetism is quite compatible 
with the Stottct (mean field) picture, as stated by sev- 
eral authorsOEJ, the T = splitting should amount to 
AE = JS. Unfortunately, the results in Figure ^| do not 
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FIG. 1: Spin dependent (solid lines up-spin, broken lines 
down-spin) bandstructure of bulk EuS calculated within a 
TB-LMTO scheme with the 4f levels treated as core states. 
The energy zero coincides with the Fermi energy. 

fully confirm this simple assumption but rather point to 
an energy-dependent exchange splitting. The indicated 
shifts of the lower edge and of the center of gravity lead 
to different J values: 

J = J(edge) = O.lfeV; J — J(center) = 0.23eV (9) 

In the following we will use both values for the a bit over- 
simplified ansatz Hdf in Eq.(Q) to compare the slightly 
different consequences. . 

It is a well-known fact (see Figure [j] irc£l and further 
references there) that the KLM can exactly be solved for 
the ferromagnetically saturated (T = 0) semiconductor. 
It is found that the f spectrum is rigidly shifted towards 
lower energies by the amount of \ JS 1 while the J, spec- 
trum is remarkably deformed by correlation effects due 
to spin exchange processes between extended 5d and lo- 
calized 4f states. We cannot switch off the interband 
exchange interaction Hdf in the LDA code, but we can 
exploit from the exact T — result that it leads in the 
T spectrum only to an unimportant rigid shift. So we 
take from the LDA calculation, which holds by defini- 
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FIG. 2: Spin-dependent density of 5d-states of EuS as func- 
tion of energy (calculated within a TB-LMTO scheme). The 
numbers are (in eV) for the lower band edges and for the 
centers of gravity. The Fermi edge is below the 5d band (see 
Fig-0). 



tion for T = 0, the \ part as the single-particle input for 
H$d in Eq.(|^). Therewith it is guaranteed that all the 
other interactions, which do not explicitly enter the KLM 
operator (|J), are implicitly taken into account by the 
LDA-rcnormalized single-particle Hamiltonian (|5|). On 
the other hand, a double counting of any decisive inter- 
action is definitely avoided. 



III. MODEL EVALUATION 

A. Magnetic Part 

Because of the empty conduction bands the magnetic 
ordering of the localized 4f moments will not directly be 
influenced by the band states. For the purely magnetic 
properties of EuS it is therefore sufficient to study exclu- 
sively the extended Heisenberg-Hamiltonian (||). While 
the exchange integrals Jy are known from spin wave anal- 
ysis (see Eq.(||)), the single-ion anisotropy constant Do 
must be considered an adjustable parameter. Via the 
magnon-Green function 



P ij (E) = ((Sf;Sr)) 



(10) 



we can calculate all desired f spin correlation functions 
by evaluating the respective equation of motion: 



EP^E) = 2h 2 (S?) % + (([S?,Hif] _ ; S. 



(11) 



Evaluation of this equation of motion requires the de- 
coupling of higher Green functions, originating from the 
Heisenberg exchange term and the anisotropy part in 



Eq.(||). For Green functions coming out of the Heisen- 
berg term we have used the so-called Tyablikow approx- 
imation, which is known to yield reasonable results in all 
temperature regions. For Green functions, which arise 
from the anisotropji_term, we use a decoupling technique 
proposed by LinescEl. Details jaf the method have been 
presented in a previous paperEH on EuO. As result one 
gets the following well-known expression for the temper- 
ature dependent local-moment magnetization: 

/cz , , (1 + 0) 2S+1 (S - 0) + ^ 2S+1 (S + 1 + 0) 
w / ~ n i5«zi / 1 ; a\2s+i [ - iZ > 



_ (1 

can be interpreted as average magnon number: 
1 



N ^ 

k 



(13) 



where E(k) is the pole of the wave vector dependent 
Fourier transform of Py (E). Some typical magnetization 
curves are plotted in Figure |§|. They differ by the value 
of the anisotropy constant Dq, which is still an undeter- 
mined parameter. When Do/ks increases from 0.01K to 
0.4K Tc rises from about 15K to 16. 9K. Regarding that 
J\, J 2 are derived from a low-temperature spin- wave fit, 
the agreement between the calculated T c s and the exper- 
imental value of 16.57KEJ is remarkably good for almost 
all applied values of Do, and best for 0.375K. Since we 
are interested above all in the electronic bulk band struc- 
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FIG. 3: 4f-magnetization as a function of temperature for var- 
ious values of the single-ion anisotropy Do . The vertical bro- 
ken line marks the used experimental value (Tc = 16.57eV). 

ture and its temperature dependence, the actual numeri- 
cal value of Do does not play the decisive role. However, 
when treating systems of lower dimensionality (films, sur- 
faces), which is planned for a forthcoming paper, then a 
finite D will be the precondition .fevr g etting a collective 
magnetic order of the spin systemr 8 rr 9 l. 

Together with ([l2]) and (|lj) practically all local spin 
correlations are derivable, as e.g.: 

(14) 



(S-S + ) = 2h{S z )<j) 
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((S z ) 2 ) = 7i 2 S(S+l)(S z )(l + 20) (15) 



((s z ) 3 ) = n 3 s(s + 1)0 + n 2 (s z ) (S(s + i) + . 

-ft ((S*) 2 ) (1 + 20) 



(16) 



These and similar terms are responsible for the temper- 
ature dependence of the electronic selfenergy. 



B. Electronic Part 

The inspection of the electronic part starts from the re- 
tarded Green function ((cj mCT ; ct m , CT )).E or its wave- vector 
dependent Fourier transform: 



G kCT (£) = 



n 



El - T k - E kCT (£) 



(17) 



Here 1 represents the (M x M) identity matrix, where 
M is the number of relevant subbands. The elements 
rpmm Q £ hopping matrix are the Fourier transforms 
of the hopping integrals in Eq. (§) , while the elements of 
the selfenergy matrix are introduced by 



[cim„H df ]_;c} mta )) E s^JT {E)GT 3 r (E) 

Im" 

(18) 

with subsequent Fourier transformation. 

To get explicitly the selfenergy elements in Eq. ( O ) we 
evaluate the commutator on the left-hand side what pro- 
duces two higher Green functions: 



rr k TAE) = {{s*c kmr7 ;c] 



rpmm 
ikja 



(E) 



(19) 



(20) 



r arises from the Ising type interaction in the d-f inter- 
action term (M) and F from the spin exchange partial 
operator {S^= S+<-): 



[c-ima, Hdf]_ i c ]m'o-// — i z a^iij<r 



(21) 

(z a = (5^1 — Sal). Exploiting already the fact that the 
EuS conduction band is empty we encounter the following 
equations of motion of the higher Green functions ( |l9| ) 
and (|2T 



E {ESkiSmm" - Ur") r^'l" 1 ' (E) = (22 

Ira" 

h(S Z )6 k j5 mm > + U S- [c k ,na,Hdf}_ |c] m 



E [EhlSmm" ~ 1^"*") ^: m ' (25) = (23) 
lm" 

($*tS~ + S^Sf) [c fcm _ CT , ; c]mv)) B 



On the right-hand side of these equations appear further 
higher Green functions which prevent a direct solution 
and require an approximative treatment. That shall be 
different for the non-diagonal terms (i ^ k) and the di- 
agonal terms (i = k), because the strong intraatomic 
correlations due to the on-site interaction (0) have to be 
handled with special care. For i ^ k a self-consistent 
selfenergy approach is appiied^i^jaida, tas been tested in 
numerous previous papers^EjaOEjEil. It simply con- 
sists in treating the commutators in ( p2| ) and p3|), re- 
spectively, in formal analogy to the definition equation 
(O) for the selfenergy: 



i [Ckma, Hdf]_ c j m 'a 



E E 

Ira" 



mm j^rn m ( t?\ 
itu 1 ilia 



(24) 



(((S^Sr + S al S+) [c km -a, H dS ]_ ■ c] m , a )) E (25) 



yimm Tpm m 
il — a iljcr 



(E) 



For the diagonal terms (i = k) a moment technique is 
used that takes the local correlations better into account. 
For this pur pos e we explicitly evaluate the commutators 
in Eqs.([22|),(|23|) obtaining then, as usual, further higher 
Green functions. In the first step these new functions 
are reduced to a minimum number by exploiting that all 
functions, arising from the I sing- equation ([22]), can rigor- 
ously be expressed by linear combinations of those which 
come out of the spinflip- equation (23). For a decoupling, 
the latter are then written as linear combinations of sim- 
pler functions that are already involved in the hierarchy 
of Green functions. The choice, which kind of simpler 
functions enter the respective ansatz, is guided by ex- 
actly solvable limiting cases (ferromagnetic saturation, 
zero-bandwidth limit, S — ^). We illustrate the proce- 
dure by a typical example: 



D%™'(E) = ((A i<7 c im<7 ;c+ m , (r )) 



A iCT =5^SrS+ + 5„tS?S, 



a I ' 



(26) 



For S = h it holds rigorously: 



D?™ (E) = -G*£ (E) - (<5 ffT - 5 ai )T%™ (E) (27) 

This relation is valid for all temperatures. On the other 
hand, in ferromagnetic saturation ((S z ) = S) the same 
function reads for all spin values: 

D^'(E) = SG^'(E) - (6 al - 8 a[ )T^(E) (28) 

Eqs. (^) and (^8|) clearly suggest the following ansatz 
for the general case: 

D%?(E) = a™ m 'G™f(E) + [3™ m 'r™'(E) (29) 
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In order to fix the coefficients a™" 1 and /3™" 1 , we now 
calculate the first two spectral moments of each of the 
three Green functions in (p9|), and that exactly and in- 
dependently from the respective Green function. The 
diagonal parts of all other functions, appearing on the 
right-hand sides of (H|) and (|23|), can be elaborated anal- 
ogously. 

By these manipulations we arrive at a closed system of 
equations for the selfenergy matrix elements X™™ (E), 
which can numerically be solved. Via the spectral mo- 
ments, used for the various ansatze such as (^9|), a set of 
local spin correlations as those in Eqs.( ^2| , |l4|J15| , [l^ ) en- 
ter the procedure. They are mainly responsible for the 
temperature-dependence of the electronic selfenergy. 

To get a first impression of correlation effects in the 
electronic structure of EuS we have evaluated our com- 
plex theory for T=0K. The Q-DOS results are exhib- 
ited in Figure H. As explained and tentatively justi- 
fied before Eq.(p|) we use two different values for the ex- 
change coupling J. The f Q-DOS is unaffected by the 
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FIG. 4: The same as in Fig. |, but in addition the T = 
results of our combined many-body / first principles theory. 
The latter has been performed for two different values of the 
exchange coupling J (dotted line J = O.lleV; broken line 
J — 0.23eV). The up-spin curves have been shifted both to 
coincide at the lower edge with the LSDA curve. That demon- 
strates that at T — correlation effects appear in the down- 
spectrum only. 

actual value of J and coincides with the respective LDA 
curve, when we compensate, as done in Figure ^, the 
unimportant rigid shift (— hJS). So our approach ful- 
fills the exact (T = 0, a =|)-limit. The slight deviations, 
seen in the upper part of Figure ^, are exclusively due 
to the numerical rounding procedure. A 'posteriori this 
fact demonstrates once more that our above-described 
method for implementing the LDA input into the many- 
body model calculation definitely circumvents the often 
discussed double counting problem. As explained in Sec- 
tion 2 we succeeded in this respect because for the special 



case of a ferromagnetically saturated semiconductor the 
t spectrum is free of correlation effects which stem from 
the interband exchange Hdf- 

The lower half of Figure [i] demonstrates that correla- 
tion effects do appear, even at T — OK, in the j spectrum. 
Besides a band narrowing, they provoke strong deforma- 
tions and shifts with respect to the LDA result. Here the 
influence of the different J values from Eq.(||) is quite 
remarkable. For getting quantitative details of the EuS- 
energy spectrum a proper choice of the exchange constant 
is obviously necessary. The lower value (J = O.lleV) is 
appropriate when we are mainly interested in the lower 
band-edge region. However, in the middle of the band, 
around the center of gravity, J = 0.23eV is surely the 
better choice. 



IV. EUS ENERGY SPECTRUM 

The main focus of our study is the temperature de- 
pendence of the 5d-encrgy spectrum of the ferromagnetic 
semiconductor EuS. The 5d bands are empty, except for 
the single test electron, so that the T-dependence must be 
exclusively caused by the exchange coupling of the band 
states to the localized magnetic 4f states. Figure |^ and 
Figure ^ display the quasiparticle densities of states for 
five different 4f magnetizations, i.e. five different temper- 
atures (Figure |). The Q-DOS in Figure || are calculated 
for J = O.lleV. One sees that the lower edge of the "f 
Q-DOS performs a shift to lower energies upon cooling 
from T — Tq down to T = OK. This explains the fa- 
mous red shift of the optical absorption edge for an elec- 
tronic 4 / — 5 dt 2 transition, first observed by Busch and 
WachteiOH. We find a shift of about 0.17eV ijsee inset 
in Figure [5]) ,very close to the experimental datafflj. This 
confirms once more that J = O.lleV is a realistic choice 
for the exchange coupling constant as long as the lower 
part of the 5d spectrum is under consideration. Note, 
that our fitting procedure for the exchange constant J 
(Fig. ||, Eq. (j^)) does not predetermine the redshift. 

Since we have taken into account for our calculation 
the full band structure of the Eu-5d conduction bands, 
the symmetry of the different Eu-5d orbitals is preserved. 
The 5d bands of bulk EuS are therefore split into ti g and 
e g subbands (Figure pf) , where the t2 g bands are substan- 
tially broader (~ 7ev) than the e g bands (~ 4eV). In a 
previous ^ uo -y °f ^ ne temperature dependent EuS-band 
structuret3 the Eu-5d complex was split into five s-like 
bands by numbering for a given k vector the states from 
m=l to m=5 according to increasing single-electron en- 
ergies e m (k). All k states with an energy e m (k) then 
form the subband m. This simplified procedure does not 
respect symmetries and neglects subband hybridization, 
i.e. interband hopping T™ m for m ^ m! ', The subband 
widths W turn out to be of order 1 — 3eV being there- 
fore distinctly narrower than those in Figure |. That has 
an important consequence. Since correlation effects scale 
with the effective(!) exchange coupling they become 
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FIG. 5: Quasiparticle densities of states of the Eu-5d bands 
of bulk EuS as a function of energy for various temperatures. 
Solid lines for up-spin, broken lines for down-spin, thick line 
for T — Tc ((S z )=0). The outermost curves belong to T = 
({S z )/S = 1). They approach each other when increasing the 
temperature. The inset shows on an enlarged scale the tem- 
perature shift of the lower up-spin edge. Exchange coupling: 
J = O.lleV. 
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FIG. 6: The same as in Fig. | but for J = 0.23eV. 



for the same J more pronounced in narrower bands. That 
is why we believe-that correlations are to a certain degree 
overestimated inlla. As a consequence of the weaker effec- 
tive coupling the appearance of polaron-like quasiparticle 
branches is less likely in the present investigation. 

In Figure || the Q-DOS is plotted for the stronger ex- 
change coupling J = 0.23eV, which should be more re- 
alistic for the middle part of the spectrum, around the 
center of gravity. The temperature-influence on the spec- 
trum is more pronounced than for the weaker coupling in 
Figure Strong deformations and shifts appear, being 
not at all rigid, i.e. far beyond the mean field picture. 
However, not surprising, the lower edge shift between 
T = Tq and T = OK comes out too strong. The calcu- 
lated red shift of 0.27eV,fubstantially exceeds the experi- 
mental value of 0.167eVH. As mentioned above, the lower 
part of the spectrum is better described with J = O.lleV. 

While the Q-DOS refers to the spin resolved, but an- 
gle averaged (inverse) photoemission experiment, the k 
dependent spectral density is the angle resolved coun- 
terpart. From Sk a {E) we derive the quasiparticle band 
structure (Q-BS). Figures ^ and || represent as density 




W L 



X WK 



FIG. 7: Spin-dependent quasiparticle bandstructure of the 
Eu-5d bands of bulk EuS for different 4f magnetizations 
(S z )/S. Exchange coupling: J = O.lleV. 



plots the spectral density for some symmetry directions. 
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T — ► T c ((5 Z ) — > 0) the vanishing 4f magnetization re- 
moves the induced spin asymmetry in the 5d subbands. 



W L 



X WK 



FIG. 8: The same as in Fig. §, but for J = 0.23eV. 



The degree of blackening is a measure of the magnitude of 
the spectral function. Figure [?] holds for J = 0.1 leV and 
Figure || for J = 0.23eV. In both situations the | spec- 
trum in case of ferromagnetic saturation coincides with 
the dispersions obtained from the LDA calculation. In 
the weak coupling case (Figure ^) the temperature influ- 
ence is mainly a shift of the total spectrum. The induced 
exchange splitting reduces with increasing T and disap- 
pears at T = Tq- Correlation effects are more clearly 
visible in the case of J = 0.23eV (Figure |J). They man- 
ifest themselves above all in lifetime effects. Great parts 
of the dispersions are washed out because of magnon ab- 
sorption (emission) of the itinerant f (J,) electron with 
simultaneous spinflip. In ferromagnetic saturation a j 
electron has no chance to absorb a magnon because there 
does not exist any magnon. Therefore the dispersion ap- 
pears sharp representing quasiparticles with infinite life- 
time. On the other hand, the J. electron has even at 
T = OK the possibility to emit a magnon becoming then 
a f electron. Therefore correlation effects are already 
at T = OK present in the j spectrum. For finite tem- 
peratures, finite demagnetizations, magnons are available 
and absorption processes provoke quasiparticle damping 
in the f spectrum, too. The overall exchange splitting 
reduces with increasing temperatures, until in the limit 




FIG. 9: Spin-dependent spectral density Ska of the Eu-5d 
states of bulk EuS as function of energy for the same 4f mag- 
netizations as in Fig. ^| Solid lines: up-spin; broken lines: 
down-spin; thick lines: T — Tc ({S z ) = 0)). Exchange cou- 
pling: J = 0.23eV; k = T. 




FIG. 10: The same as in Fig. 9 but k = L. 




FIG. 11: The same as in Fig. g but k = W. 

For a better insight into the temperature-behavior we 
have plotted in Figs. p|-p"2| for four k-points from the 
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FIG. 12: The same as in Fig. § but k = X. 



Brillouin zone (r, L, W, X) the energy dependence of 
the spectral density, and that for the same three tem- 
peratures as in Fig. ||. For the T-point we expect two 
structures according to the twofold (e g ) and threefold 
(t2 g ) degenerate dispersions. As can be seen in Fig. ^ 
well defined quasi particle peaks appear with additional 
spin split below To The exchange splitting, introduced 
via the interband coupling to the magnetically active 4f 
system, collapses for T — > Tq {Stoner-like behavior). Ob- 
viously, quasiparticle damping increases with increasing 
temperature. Similar statements hold for the spectral 
density at the L-point. In accordance with the quasipar- 
ticle bandstructure in Fig. three structures appear, the 
upper two being twofold degenerate (Fig. HO). Interest- 



ing features can be observed at the W-point (Fig. 11). 
At T = four sharp peaks show up in the f -spectrum, 
and, though already strongly damped, the same peak- 
sequence comes out in the |-spectrum. The exchange 
splitting amounts to about 0.8-0.9eV. With increasing 
temperature damping leads to a strong overlap of the 
two upper peaks, which are no longer distinguishable. 

At the X-point the spin resolved spectral density ex- 
hibits four clearly separated structures, where the upper 
belongs to a twofold degenerate dispersion (Fig. |l^). In 
the two middle structures, at least the J, -peaks are so 
strongly damped that they certainly will not be observ- 
able in an inverse photoemission experiment. Altogether, 
the 5d-spectral densities of the ferromagnetic semicon- 
ductor EuS exhibit drastic temperature-dependencies, 
what concerns the positions and the widths of the quasi- 
particle peaks. 



V. CONCLUSIONS 

We presented a method of calculating the temper- 
ature dependent bandstructure for the ferromagnetic 
semiconductor EuS. The essential point is the combina- 
tion of a many body evaluation of a proper theoretical 
model with a "first principles" band structure calcula- 
tion. The model of choice is the ferromagnetic Kondo- 
lattice model, which describes the exchange interaction 



between localized magnetic moments and itinerant con- 
duction electrons. For a realistic treatment of EuS we 
have extended the conventional KLM to a multiband ver- 
sion to account for orbital symmetry. Intraband- and 
interband-hopping integrals have been extracted from a 
tight-binding linear muffin-thin orbital band structure 
calculation to incorporate besides the normal single- 
particle energies the influence of all those interactions 
which are not directly covered by the KLM-Hamiltonian. 
By exploiting an exactly solvable limiting case of the 
KLM this combination of first-principles and model- 
calculations could be done under strict avoidance of a 
double-counting of any relevant interaction. 

The many-body part of the procedure was performed 
within a moment-conserving interpolation method that 
reproduces exactly important rigorous limiting cases of 
the model. The resulting electronic selfenergy carries a 
distinct temperature-dependence, which is mainly due to 
local 4f spin correlations. Since the band is empty, the 
KLM reduces to a simple Heisenberg model as long as 
the purely magnetic 4f properties, as e.g. the mentioned 
4f spin correlations, are concerned. The result is a closed 
system of equations which can be solved numerically for 
all quantities we are interested in. 

We have demonstrated the temperature-dependence of 
the energy spectrum of the ferromagnetic semiconductor 
EuS in terms of the 5d spectral density and 5d quasipar- 
ticle density of states. Peak positions and peak widths 
determine energies and life-times of quasiparticles, which 
have been gathered in special quasiparticle band struc- 
tures. A striking temperature-dependence of the empty 
5d-bands is observed which is induced by the magnetic 
4f-state. A well-known experimental confirmation of the 
T-dependence is the "red shift" of the optical absorption 
edge, uniquely related to the shift of the lower 5d-edge for 
decreasing temperature from T = Tq down to T = OK. 
The induced exchange splitting at T = collapses for 
T — > Tq Stoner-like, but with distinct changes in the 
quasiparticle damping (lifetime) . All these temperature- 
dependent band effects should be observable by use of 
inverse photoemission. 

We expect further insight into interesting physics by a 
forthcoming extension of our method to finite band oc- 
cupations (Gd, Gd-films). The respective single-band 
model . version has already been presented in previous 



papeaM-M— I In particular a modified RKKY has been 



usedH3 for a selfconsistent inclusion of the magnetic prop- 
erties of the not directly coupled local moments. A fur- 
ther implementation of disorder in the local moment sys- 
tem will allow to treat diluted magnetic semiconductors 
such as Gai-^Mn^As and therefore contribute to the hot 
topic spintronics. 
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